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Abstract 

The properties of ripplonic polarons in a multielectron bubble in liquid helium 
are investigated on the basis of a path- integral variational method. We find 
that the two-dimensional electron gas can form deep dimples in the helium 
surface, or ripplopolarons, to solidify as a Wigner crystal. We derive the 
experimental conditions of temperature, pressure and number of electrons in 
the bubble for this phase to be realized. This predicted state is distinct from 
the usual Wigner lattice of electrons, in that it melts by the dissociation of the 
ripplopolarons, when the electrons shed their localizing dimple as the pressure 
on the multielectron bubble drops below a critical value. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The two-dimensional (2D) electron system formed on the surface of liquid helium has 
been widely investigated, especially with regard to the formation and melting of a Wigner 
lattice [1]. An electric field pressing the electrons against the hehum surface results in 
an interaction between the electrons and the quantized modes of oscillation of the helium 
surface, the ripplons [2]. In this paper, we investigate the effects of the clcctron-ripplon 
coupling in a multielectron bubble in liquid helium and highlight the differences with the 
case of electrons on a fiat hehum surface. Such multielectron bubbles (MEBs) are spherical 
cavities in liquid hehum. MEBs are typically 0.1 //m - 100 //m in radius for 10^—10^ electrons 
[3]. The bubble's diameter is determined by the balance between the surface tension of liquid 
helium and the Coulomb repulsion of the electrons [4]. 

The electrons are not distributed throughout the volume of the bubble, but form a 
nanometer thin and effectively 2D layer anchored to the inside of the bubble surface [5]. 
Although the flat 2D electron gas has been studied extensively and new quantum states 
such as the fractional quantum Hall regime were revealed, much less effort has gone to the 
study of the spherically curved 2D electron gas such as the one present in MEBs. Spherical 
shells of charged particles appear also in a variety of other physical systems, such as fullcrcncs 
[6], metallic nanoshells [7], and charged macroscopic droplets. The main goal of this paper 
is to show that the enhanced electron-ripplon coupling in the bubble leads to a new solid 
phase, a lattice of ripplonic polarons, that is distinct from the electron Wigner lattice, and 
to investigate the properties of this phase. 

Using pressure, the surface of the MEB can be compressed to achieve 2D electron densities 
as high as 10^^ cm~^ [8], whereas flat surfaces are hmited to 2 x 10^ cm~^ due to an instabihty 
[9]. This instability is not present in multielectron bubbles [10]. As a result very large electric 
flelds exist on an electron, normal to the surface (due to all other electrons in the MEB), 
whereas for a flat surface the maximum field is around 3 kV/cm. Some typical values for 
physical variables related to the bubble are given in table I. 

The Hamiltonian of a single electron on a fiat helium surface is given by 



where p is the electron momentum operator, mg is the electron mass, and u!{q) — q 



is the ripplon dispersion relation with a ~ 3.6 x 10~^ J/m^ the surface tension of helium 
and p — 145 kg/m^ the mass density of helium. In this Hamiltonian, we restrict ourselves 

to 2D position and momentum operators, assuming that the part of the wave function 
of the electrons relating to the direction perpendicular to the surface can be factored out 
exactly. The second-quantization operators d^,aq create/annihilate a ripplon with planar 
wave number q. The electron-ripplon coupling amplitude is given by 



where E is the electric field perpendicular to the surface (the so-called 'pressing field'), and 
e is the electron charge. The pressing field pushes the electrons with a force eE towards 
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the helium surface, that acts hke a sheet with surface tension a. Note that there is a 1 cV 
barrier preventing single electrons from penetrating the helium surface. The self- induced 
trapping potential of the electron on the helium surface is manifested by the appearance 
of a dimple in the helium surface underneath the electron, much like the deformation of a 
rubber sheet when a person is pulled down on it by a gravitational force. The resulting 
quasiparticle consists of the electron together with its dimple and can be called a ripplonic 
polaron or ripplopolaron [11]. 

Hamiltonian (1) for the ripplopolarons is very similar to the Frohlich Hamiltonian de- 
scribing polarons [12]; the role of the phonons is now played by the ripplons. Methods 
suitable for the study of single polarons have been used to analyse the single ripplopolaron 
on a fiat surface [13,14]. Recently, Fratini and Qucmerais [15] have proposed a path integral 
treatment for a Wigner lattice of polarons. One of the goals of the present paper is to adapt 
their method so that it becomes suitable for the treatment of a lattice of ripplopolarons. 

In Sec. II we introduce the Hamiltonian describing a ripplopolaron in a Wigner lattice 
of ripplopolarons in a multielectron bubble. The temperature zero solution for the ground 
state of this Hamiltonian is derived in the strong-coupling case in Sec. III. The results for 
arbitrary interaction strength and arbitrary finite temperature, obtained with the variational 
path-integral technique, are presented in Sec. IV. These finite-temperature results allow us 
to investigate the melting of the ripplopolaron Wigner lattice and determine the phase 
diagram of this state in Sec. V. The results are discussed and compared with the case of 
electrons on a flat hehum surface in Sec. VI. 



II. HAMILTONIAN FOR A RIPPLOPOLARON IN A WIGNER LATTICE 

In their treatment of the electron Wigner lattice embedded in a polarizable medium such 
as a semiconductors or an ionic solid, Fratini and Quemerais [15] described the effect of the 
electrons on a particular electron through a mean-field lattice potential. The (classical) 
lattice potential Viat is obtained by approximating all the electrons acting on one particular 
electron by a homogenous charge density in which a hole is punched out; this hole is centered 
in the lattice point of the particular electron under investigation and has a radius given by 
the lattice distance d. Thus, in their approach, the anisotropy effects are neglected. A 
second assumption implicit in this approach is that the effects of exchange are neglected. 
This can be justified by noting that for the electrons to form a Wigner lattice it is required 
that their wave function be localized to within a fraction of the lattice parameter as follows 
from the Lindemann criterion [16]. As can be read from table I, the typical distance between 
electrons (the lattice parameter) is 10-100 nm. 

Within this particular mean-field approximation, the lattice potential can be calculated 
from classical electrostatics and we find that for a 2D electron gas it can be expressed in 
terms of the elliptic functions of first and second kind, E {x) and K {x\ 
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Here, r is the position vector measured from the lattice position. We can expand this 
potential around the origin to find the small- amplitude oscillation frequency of the electron 
lattice: 

hm (r) = - ^ + \m^u;iy + O {r ") , (4) 
with the confinement frequency 



<^lat — 




Although this mean-field approach may seem crude, it has the distinct advantage that the 
'phonon' frequency ujiat of the electron lattice corresponds closely to the longitudinal plasmon 
frequency that can be derived using an entirely different approach based on a more rigorous 
study of the modes of oscillations of both the bubble and the charge distribution on the 
bubble surface. This frequency lies typically in the THz range and the lattice parameter d 
in MEBs ranges roughly from 10 to 100 nm. From this, and from the succesful application of 
this mean-field approach to polaron crystals in solids, we conclude that the approach based 
on that of Fratini and Quemerais describes the infiuence of the other electrons well in the 
framework of small amplitude oscillations of the electrons around their lattice point. The 
(modified) Lindemann melting criterion suggests that the lattice will melt when the electrons 
are on average displaced more than ca. 10% from their lattice position; thus in the regime 
of interest the Pratini-Quemerais approach is applicable. In the mean-field approximation, 
the Hamiltonian for a ripplopolaron in a lattice on a locally flat helium surface is given by 

H^£^ + Viat{r) + Yl hw{q)dXd^ 

+ Y^M,e-^^-^(a^ + dt^, (6) 
q 

where f is the electron position operator. 

Now that the lattice potential has been introduced, we can move on and include effects of 
the bubble geometry. If we restrict our treatment to the case of large bubbles (with > 10^ 
electrons) such as those already experimentally observed [3], then both the ripplopolaron 
radius and the inter-electron distance d are much smaller than the radius of the bubble 
R},. This gives us ground to use the locally flat approximation using the auxiliary model 
of a ripplonic polaron in a planar system described by (6), but with a modified ripplon 
dispersion relation and an modified pressing field. In appendix A we provide a more detailed 
description of the electron-ripplon interaction on a sphere and discuss how the locally flat 
approximation (6) is linked to the exact expressions for the curved surface. Essentially, we 
flnd for the modifled ripplon dispersion relation in the MEB: 

V P pRh 

where i?b is the equilibrium bubble radius which depends on the pressure and the number of 
electrons [8]. The bubble radius is found by balancing the surface tension and the pressure 
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with the Coulomb repulsion [8]. At zero pressure, it scales as N"^/^ and in the pressure- 
dominated regime {p ^ cr/Rh) it scales with pressure as p~^/^. Some typical values are 
given in table I. The modified electron-ripplon interaction amplitude in an MEB is given by 



The effective electric pressing field pushing the electrons against the helium surface and 
determining the strength of the electron-ripplon interaction is 

Ne 

E = -^e.. (9) 

Some typical values for the pressing field and the dimensionless coupling constants are given 
in Table I. Note that the electric field in the bubble is larger than the typical pressing fields 
(of the order of 10^ — 10^ V/cm) apphed on electrons on a flat hehum surface. Thus, the 
electron-ripplon coupling will be stronger in the multielectron bubble. The modified ripplon 
dispersion relation (and the dependence of the bubble radius on the number of electrons and 
the pressure) was studied in more detail in Ref. [8], and the stability of the multielectron 
bubble against surface deformational modes was investigated in detail in Refs. [5,10]. These 
studies concluded that even though a large effective electric pressing field is present at zero 
pressure, the bubbles can be stable in contrast to the flat surface which can only sustain a 
moderate electric pressing field. 

The crucial differences that exist between the case of a ripplopolaron in the multielectron 
bubble and on the ffat surface (and that are preserved in the locally ffat approximation) are 
(i) the electric pressing field E is stronger than that typically realised for electrons on helium 
films (see table I) and thus the electron-ripplon coupling is enhanced as compared to the 
normal film; (ii) the interaction energy arising from the change in polarisation of the helium 
due to the displacement of the electron has a similar form, but is much weaker and can be 
neglected. In addition the electric field, and thus the electron-ripplon coupling increases as 
the bubble radius is decreased. Thus, pressurizing the bubbles, which decreases the radius, 
also increases the electron-ripplon coupling strength (rougly as R^^)- In the high- pressure 
regime (p ^ a/Rb), the bubble radius goes as p~^^'^ and thus the electron-ripplon coupling 
increases as y^. The pressure provides a 'tuning knob' to set the electron-ripplon interaction 
strength at a desired level. 

III. THE RIPPLOPOLARON WIGNER LATTICE AT TEMPERATURE ZERO 



To gain insight into the nature of the Wigner solid of ripplopolarons, we will analyse 
Hamiltonian (6) first in the strong- coupling approach. In the next section, the more general 
and more accurate Feynman variational path-integral method will be applied, generalizing 
the results of this section to finite temperature. Noting that the frequency associated with 
the electron's motion, ooiat, is typically several orders of magniture larger than the frcqiiency 
associated with the ripplons [17], a;(g), we can safely make the product ansatz for the wave 




(8) 



A. Ground state 
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function of the ripplopolaron in the lattice: |^) = \ipe) Here |0) is the factor of the 
wave function that contains the ripplon coordinates, and is the electronic part of the 
wave function. For small-amplitude oscillations of the electrons around their lattice site, 
the lattice potential Viat is well approximated by a parabolic potential (4), so we choose a 
Gaussian trial wave function for the electronic part: 



7rV2a 



(10) 



In this trial wave function, the variational parameter is a, the width of the electron wave 
function. Taking the expectation value of Hamiltonian ( 6) with respect to this electronic 
part of the wave function yields: 



+ 



(11) 



The ripplonic part of (^/^e|-^|^e) represents a displaced harmonic oscillator and can be 
rewritten as 



+ 



E 



|M,|2e-"'^'/2 



fiuj(q) 



hu!{q) 



fiu>{q) 



(12) 



The ground state of the displaced (2D) harmonic oscillator at temperature zero has energy 
hu!{q), independent of the variational parameter a. To find the variational optimal value of 
a, we minimize the ripplopolaron energy 



E 



q ^^{q) 

The sum over momenta can be rewritten as an integral, remembering 

f (Pq 
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The lower limit appears since the largest wavelength available is l/Rb- We checked that 
the final results do not depend crucially on the value of this naturally occurring cut-off. A 
dimensionless integration variable q = qRj, is introduced. The ground state energy of the 
ripplopolaron can then be evaluated analytically: 
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where F is the incomplete gamma function. Fig. 1 shows the result for the variational 
parameter a as a function of number of electrons and pressure in the multielectron bubble. 
If this figure, a is expressed relative to the interelectron distance ^AnRl/N. 



B. Dimple shape 

The ripplonic part of the Hamiltonian (12) represents oscillations of the helium surface no 
longer around the original bubble surface, but around a new, displaced equilibrium surface. 
This displacement of the helium surface is the dimpling. Underneath each electron, a dimple 
appears. The new equilibrium surface, described by the function u{r) (cf. appendix, Eq. 
Al), can be found by using the canonical relation between the surface displacement operator 
and the ripplon creation and annihilation operators: 

= + ^^^^ 

and evaluating 

«(r) = ^(^|gq|^)e^'»-. (18) 
q 

The result is given by 

H 

a 

In the limiting case of a large bubble, this result corresponds to that of Shikin and Monarkha 
[11] for electrons on a flat helium surface; the role of the capillary constant is played by 
p/{aRb). Fig. 2 shows, for a bubble with N — 10^ electrons, at different pressures the 
shape of the dimpled surface. Several dimples are shown - above the center of each dimple 
an electron is present. The dotted curve represents the undimpled u{r) — surface; the 
curvature of the bubble surface is visible in this curve. The electrons are separated by the 
interelectron distance d = ^JAttRI/N. As the pressure increases, the radius of the bubble 
decreases. Since the number of electrons is fixed the electric pressing field increases, making 
on its turn the electron-ripplon coupling larger. This results in deeper, narrower dimples. 
Note that while the deformation here can be several angstroms, for a flat surface on bulk 
helium the maximum deformation of a dimple is less than one angstrom [18]. Also for 
electrons on a thin helium film above a dielectric substrate, the dimple depth can reach 
several angstrom [19]. 
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IV. THE RIPPLOPOLARON WIGNER LATTICE AT FINITE TEMPERATURE 



The simple but intuitive approach of the previous section describes the system in the 
hmit of temperature zero. To study the ripplopolaron Wigner lattice at finite temperature 
(and for any value of the electron- ripplon coupling), we use the variational path-integral 
approach [20]. This variational principle distinguishes itself from Rayleigh-Ritz variation in 
that it uses a trial action functional Striai instead of a trial wave function. 

The action functional of the system described by Hamiltonian (6), becomes, after elimi- 
nation of the ripplon degrees of freedom. 





Jdrj daG^^,){T - a)e*'i-['-(^)-'-('^)], (20) 



X 





with 



cosh[z/(|r-a|-^/?/2)] 
^^^^ - ^) = sinh(/?W2) ■ 

In preparation of its customary use in the Jensen-Feynman inequality, the action functional 
(20) is written in imaginary time t = ir with f3 = 1/(/cbT) where Tis the temperature. 
Following an approach analogous to that of Fratini and Quemerais for a lattice of polarons 
in a semiconductor [15], and to that of Devreese et al. for N polarons in a quantum dot 
[21], we introduce a quadratic trial action of the form 



Striat ~ ^ 







' dr I daG^{T - (j)r(r) • v{a). (22) 



4:% 



where M, w, and Q are the variationally adjustable parameters. This trial action corresponds 
to the Lagrangian 

from which the degrees of freedom associated with R have been integrated out. This La- 
grangian can be interpreted as describing an electron with mass rUe at position r, coupled 
through a spring with spring constant k to its lattice site, and to which a fictitious mass M 
at position R has been attached with another spring, with spring constant K. The relation 
between the spring constants in (23) and the variational parameters w, Q is given by 

w = ^jK/m,, (24) 
n = ^J{K + K)/me. (25) 
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Based on the trial action Striai, Feynman's variational method allows one to obtain an 
upper bound for the free energy F of the system (at temperature T) described by the action 
functional S by minimizing the following function: 



trial I 



(26) 



with respect to the variational parameters of the trial action. In this expression, Fq is the 
free energy of the trial system characterized by the Lagrangian £o) 13 = l/{khT) is the 
inverse temperature, and the expectation value {S — Striai) is to be taken with respect to 
the ground state of this trial system. The evaluation of expression (26) is straightforward 
though lengthy. We find 



In 



2 sinh 



13 



\/7re 
D 
1 



2 sinh 
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+ — In 2 sinh 

h 2 / 

- afQi coth \^ 



-(f/{2D) 



dqq\Mg 



' / dr 
Jo 



2D^ 

cos\i[uj{q){T -h(3/2)] 
smh.[f3huj{q)/2] 



(27) 



X exp 



2me ^ 
.7=1 



2 cosh(ftOj 13/2) -cosh [nVLj (r-/3/2)] 

3 Qj smh{hCljf3/2) 



In this expression, /q and Ii are Bessel functions of imaginary argument, and 

h 2 ^ 



^ = — E ^ coth {hnj(3/2) , 

Trie -^^ \Lj 



(28) 



ai 



w 



w 



(29) 



Finally, fti and 0,2 are the eigenfrequencies of the trial system, given by 



=0 + ± ^(02 _ ^2)2 ^ 4^K/{Mme 



(30) 



Optimal values of the variational parameters arc determined by the numerical minimization 
of the variational functional F as given by expression (27). As the reader may notice, the 
result of the variational path-integral method is slightly less intuitive than that of the strong- 
coupling approach of the previous section, nevertheless it is much more general and will allow 
us to introduce temperature to examine the melting of the Wigner lattice of ripplopolarons 
in the next section. 
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V. MELTING OF THE RIPPLOPOLARON WIGNER LATTICE 



The Lindemann melting criterion [16] states in general that a crystal lattice of objects (be 
it atoms, molecules, electrons, or ripplopolarons) will melt when the average motion of the 
objects around their lattice site is larger than a critical fraction So of the lattice parameter d. 
It would be a strenuous task to calculate from first principles the exact value of the critical 
fraction 60, but for the particular case of electrons on a helium surface, we can make use 
of an experimental determination. Grimes and Adams [22] found that the Wigner lattice 
melts when F = 137 ± 15, where F is the ratio of potential energy to the kinetic energy 
per electron. In their experiment, the electron density varied from 10^ cm~^ to 3 x 10^ 
cm~^ while the melting temperature T,. varied from 0.23 K to 0.66 K. At temperature T the 
average kinetic energy in a lattice potential Viat is 



and the average distance that an electron moves out of the lattice site is determined by 



Prom this we find that for the melting transition in Grimes and Adams' experiment [22], the 
critical fraction equals Sq ~ 0.13. This estimate is in agreement with previous (empirical) 
estimates yielding Sq ~ 0.1 [23], and we shall use it in the rest of this paper. 

The unmodified Lindemann criterion as stated above cannot be applied to an infinite 
layer of electrons on helium at non-zero temperature, because (when a thermal occupation 
of the ripplon modes is assumed) a straightforward calculation of the average distance that 
an electron moves out of its lattice site yields a divergent result. This divergence is closely 
related to Hohenberg's theorem forbidding Bose-Einstein condensation in 2D. Therefore, 
many authors rely on a modified Lindemann criterion [24] that considers the average dis- 
tance between two nearest neighbors instead of the average distance of a lattice resident 
from its lattice site. However, for the current geometry this modification is unnecessary: the 
multielectron bubble is a finite and confined system, for which considerations based on Ho- 
henberg's theorem do not apply. Hence, we shall use the unmodified Lindemann criterion to 
study the melting of the ripplopolaron lattice. In practice, we see that the above-mentioned 
divergence is not present because there is a natural cut-off wavelength for the ripplons: the 
lowest ripplon mode on a sphere corresponds to an ^ = 1 spherical harmonic, to which a 
characteristic wavelength of the order of 1/Rb can be associated. We have checked that the 
results do not depend on the precise value of the cut-off wavelength X/Rb with A on the 
order of 1. 

Within the approach of Pratini and Quemerais [15], the Wigner lattice of (ripplo)polarons 
melts when at least one of the two following Lindemann criteria are met: 




(31) 




(32) 




(33) 




(34) 
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where p and Rcms are, respectively, the relative coordinate and the center of mass coordinate 
of the model system (23): if r is the electron coordinate and R is the position coordinate of 
the fictitious ripplon mass M, this is 

mer + MR 

Rcms = rrir'^ p = r - R. (35 

The appearance of two Lindemann criteria takes into account the composite nature of (rip- 
plo)polarons. As follows from the physical sense of the coordinates p and Rcms, the first 
criterion (33) is related to the melting of the ripplopolaron Wigner lattice towards a rip- 
plopolaron liquid, where the ripplopolarons move as a whole, the electron together with 
its dimple. The second criterion (34) is related to the dissociation of ripplopolarons: the 
electrons shed their dimple. 

The path-integral variational formalism outlined in the previous section allows to calcu- 
late the expectation values (Rc^^) and (p^) with respect to the ground state of the varia- 
tionally optimal model system. We find 

/■R2 \ ^ ^ 

\-^cm./ [^2(J^2 ^ _ (q2 _ q2) 

X [q^(Q? - w"^) coth(;iQi/3/2)/Qi 

+n\{w'^ - nl) cot\i{hVt2l3/2)/Vt^ , (36) 

n 



me - ni) (fi2 _ yj2^ (yj2 _ f]2) 



X 



+Q^(Qf - w^) coth(;iQ2/3/2)J . (37) 

The procedure to find whether the Lindeman criteria are fulfilled is then as follows: first the 
optimal values of the variational parameters are obtained by minimization of the free energy 
(27), and then these optimal values are substituted in (36), (37). Numerical calculation 
shows that for ripplopolarons in an MEB the inequality f2i ^ w is fulfilled {w/Qi ~ 10~^ 
to 10~^) so that the strong-coupling regime is realized, in agreement with the results of Sec. 
II. Owing to this inequality, we find from Eqs. (36), (37) that 

^2 



Kms)«{p')- (38) 

So, the destruction of the ripplopolaron Wigner lattice in an MEB occurs through the 
dissociation of ripplopolarons, since the second criterion (34) will be fulfilled before the first 
(33). The results for the melting of the ripplopolaron Wigner lattice are summarized in the 
phase diagram shown in Fig. 3. For every value of A^, pressure p and temperature T in an 
experimentally accessible range, this figure shows whether the ripplopolaron Wigner lattice 
is present (points above the surface) or molten (points below the surface). Below a critical 
pressure (on the order of 10^ Pa) the ripplopolaron solid will melt into an electron liquid. 
This critical pressure is nearly independent of the number of electrons (except for the smallest 
bubbles) and is weakly temperature dependent, up to the helium critical temperature 5.2 
K. This can be understood since the typical lattice potential well in which the ripplopolaron 
resides has frequencies of the order of THz or larger, which correspond to ~ 10 K. 
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VI. DISCUSSION 



In the previous section wc have estabhshed that the ripplopolaron Wigner lattice will not 
melt into a liquid of ripplopolarons, but rather melt through dissociation of the composite 
quasiparticle that is the ripplopolaron. The absence of a ripplopolaron liquid phase can be 
understood intuitively from the fact that the ripplon frequencies (typically GHz) are several 
orders of magnitude smaller than the electron frequencies in the lattice potential (typically 
THz). In order to create a liquid of ripplopolarons, the ripplopolarons have to move with an 
average velocity large enough to keep the ripplopolaron lattice molten. This motion has to 
be of the entire object, namely the electron and its dimple. But, at the velocities required 
to keep the ripplopolaron liquid from freezing into a lattice, the dimples cannot follow the 
electrons. Thus, ripplopolarons only exist in a crystallized state. 

The present treatment does not allow us to derive the stucture of this lattice - the 
mean-field approximation made for the lattice potential prohibits this. The problem of the 
exact lattice structure is comphcated by the topology of the surface [25]: unhke for a flat 
surface, it is impossible to tile a sphere with a triangular lattice; frustration of the lattice in 
the form of point defects is unavoidable, providing nucleation points for melting the lattice. 
The problem of placing classical point charges on a sphere was first considered by Thomson 
[26] and was recently reconsidered for localized electrons in multielectron bubbles [27]. 

The present treatment does allow us to study also the electron Wigner lattice, by putting 
e|E| = in the above results, thus switching the electron-ripplon coupling off. A Wigner 
lattice of electrons is to be distinguished from a Wigner lattice of ripplopolarons. The 
lattice of ripplopolarons on the one hand melts through dissociation of the ripplopolarons, 
and this melting line is almost temperature independent. The lattice of electrons on the 
other hand melts though either classical thermal motion (when the temperature reaches 
a melting temperature of about 0.5 K), or through quantum melting when the density of 
electrons is large enough so that the extent of the zero-point motion becomes comparable 
to the lattice parameter. In the Wigner lattice of the ripplopolarons on the one hand, the 
particles are localized by the self-induced polaronic trapping potential (the dimple) due to 
the electron-ripplon interaction. In the Wigner lattice of electrons, the electrons are localized 
through the Coulomb interaction between the electrons. Finally, the region in phase space 
where the ripplopolaron Wigner lattice resides is different from the region where the electron 
Wigner lattice is found. This is illustrated in Fig. 4, where the phase diagram drawn in 
Fig. 3 is extended to huge bubbles (approaching the flat surface geometry). In the corner of 
largest N (N > 10^, Rb^l nim, Ug ^ 10^ — 10^° cm~^) and lowest pressure (p < 0.1 Pa), we 
find that an electron Wigner lattice (without individual dimples) can still be formed below 
T = 0.4 K. Thus, the electron Wigner lattice is recovered and the melting temperature 
derived from our treatment is in agreement with the experimentally observed temperature 
[22]. Our calculations show that, as the bubble is compressed, the electron Wigner lattice will 
quantum melt because of the increased density of electrons. The region of phase space where 
the electron Wigner lattice is present is separated from the region where the ripplopolaron 
Wigner lattice is present by a region where the predicted phase is an electron liquid. The 
ripplopolaron liquid phase, as mentioned before, does not exist. Since, as mentioned in the 
previous paragraph, our method does not allow us to study the crystal structure, we cannot, 
in the electron Wigner lattice phase, distinguish between the crystalline and the hexatic 
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phase [28]. 

The new phase that we predict, the ripplopolaron Wigncr lattice, will not be present 
for electrons on a flat helium surface. At the values of the pressing field necessary to 
obtain a strong enough electron-ripplon coupling, the flat helium surface is no longer stable 
against long- wavelength deformations [9]. Multielectron bubbles, with their different ripplon 
dispersion and the presence of stabilizing factors such as the energy barrier against fissioning 
[10], allow for much larger electric fields pressing the electrons against the helium surface. 
The regime of N,p,T parameters suitable for the creation of a ripplopolaron Wigner lattice 
lies within the regime that would be achievable in recently proposed experiments aimed 
at stabihzing multielectron bubbles [29]. The ripplopolaron Wigner lattice and its melting 
transition might be detected by spectroscopic techniques [22,30] probing for example the 
transverse phonon modes of the lattice [31]. 

VII. CONCLUSIONS 

In this paper, we investigate the properties of ripplopolarons in a multielectron bubble 
in helium using path-integral methods similar to those developed for a lattice of polarons 
[15]. Expressions are derived for the free energy and the dimple shape of ripplopolarons in 
a Wigner lattice in a multielectron bubble, as a function of temperature, externally applied 
pressure and number of electrons in the bubble. We find that, owing to the difference in 
the ripplon and longitudinal plasmon frequencies [17], the ripplopolarons exist only in a 
Wigner crystallized state. This state differs from the Wigner lattice of electrons, in that the 
electrons in the ripplopolaron Wigner lattice are localized by the electron-ripplon interaction 
rather than the Coulomb repulsion, and in that the melting occurs through the dissociation 
of the ripplopolarons. As electron-ripplon interaction is weakened (for example by reducing 
the externally applied pressure on the multielectron bubble) the electrons can shed their 
localized dimple and the ripplopolaron Wigner state is destroyed. The melting transition is 
shown to occur in a region of phase space that is accessible to recently proposed experiments 
for stabilizing multielectron bubbles. 
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APPENDIX A: ELECTRON-RIPPLON INTERACTION IN MEBS 

1. Ripplon dispersion on a sphere 

The modes of oscillation of the bubble surface will be called 'ripplons' in analogy with 
the surface oscillations on a fiat helium surface. In general, the deformed bubble surface can 
be described by a function R{Q) which gives the distance, from the center of the bubble, 
of the bubble surface in the direction determined by the spherical angles fl = {6,(j)}. The 
deformation u{fl) from spherical symmetry can be expanded in spherical harmonics Y^ot(O): 
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^Rb + ^Rb + Yl Qe,mYi,mi^). (Al) 

In this expression, is the angle-averaged radius, and Qim is the amphtude of the de- 
formational mode corresponding to the spherical harmonic Y^^m- In the summation, we 
abbreviate J2tm = S^^i I]m=-£- The equilibrium bubble radius is found by balancing the 
surface tension and pressure terms with the Coulomb repulsion: Rb satisfies 

2.R, + pRl = ^^, (A2) 

with (T fa 3.6 X 10~^ J/m^ the surface tension of helium, e — 1.0572 the dielectric constant 
of helium, p the difference in pressure outside and inside the bubble, e the electron charge 
and N the number of electrons in the bubble. Expanding the energy of the bubble up to 
second order in allows to derive the frequency of oscillation of a particular mode of 
deformation [8]: 




a(£^ + £ + l)+pi?,-£^|^ 



(A3) 



where p 145 kg/m^ the mass density of hehum. Thus, for small amphtude deformations 

of the bubble, we find that the shape of the bubble oscillates with frequencies given by (A3). 
Taking the bare ripplon frequency (without the effect of the interaction with electrons) and 
putting Rb ^ oo (with £/Rb = q a. constant) we find 




<ZM = \l -Q' + ^Q- (A4) 



This dispersion relation corresponds to the ripplon dispersion on the flat surface, with the 
difference that there is in our dispersion relation no gravitational term, but a term related 
to the pressure on the bubble. 



2. Electron ripplon interaction in the MEB 

The interaction energy between the ripplons and the electrons in the multielectron bubble 
can be derived from the following considerations: (i) the distance between the layer electrons 
and the helium surface is fixed (the electrons find themselves confined to an effectively 2D 
surface anchored to the helium surface [5]) and (ii) the electrons are subjected to a force 
field, arising fromt the electric field of the other electrons. For a spherical bubble, this 
electric field lies along the radial direction and equals 



Ne 



E=-^er. (A5) 



A bubble shape oscillation will displace the layer of electrons anchored to the surface. The 
interaction energy which arises from this, equals the displacement of the electrons times the 
force eE acting on them. Thus, we get for the interaction Hamiltonian 
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Here u{Vt) is the radial displacement of the surface in the direction given by the spherical 
angle VL; and (ij is the (angular) position operator for electron j. The displacement can be 
rewritten using (Al) and we find 

Hir^t = E e|E| E QimY^mi^j). (A7) 
3 i,m 

Using the relation 



Qe,m = (-Ij^-^-l-^l^/yH^i^l^Am + (A8) 
the interaction Hamiltonian can be written in the suggestive form 

Hint = E E Mi^rnYi,m{^j){di,m + a^-J , (A9) 

£,Tn j 

with the electron-ripplon coupling amplitude for a MEB given by 



Me,m = (-1) 



{m-\m\)/2-_ 



h{i + i) 

2Rl\ 2pRlue 



(AlO) 



3. Locally flat approximation 

Substituting M^^^ into (A9), we get 



int 



EE 

t,m j 



X 



^(m-|m|)/2 



Ne^ hiC + l) 

yg,m(^j) 

Rh 



(All) 



+ a 



l,—m) 



In this expression, we consider the limit of a bubble so large that the surface becomes fiat 
on all length scales of interest. Hence we let i?^ — > oo but keep t/Rb = q a constant. This 
means we have to let ^ — > oo as well. In this limit, 



(■^°° TTV sm 6 



sin[(£ + l/2)e + 7r/4]. 



(A12) 



and Y^^^id) varies locally as a plane wave with wave vector q = i/Rb- The wave function 
Yi^m[flj)/Rb is furthermore normalised with respect to integration over the surface (with 
total area AnRl). Thus, we get in the locally flat approximation 



Hint — EE 



hq 



3«q.rj 



2Rl V 2fxv{q) 



(A13) 
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or 



This corresponds in the hmit of large bubbles to the interaction Hamiltonian expected for a 
flat surface. 



16 



REFERENCES 



* Also at: Technische Universiteit Eindhoven, Eindhoven, The Netherlands. 
[1] Two-dimensional electron systems on helium and other cryogenic substrates (ed. E.Y. 

Andrei, Kluwer Acad. Publ., Dordrecht, The Netherlands, 1997). 
[2] V.B. Shikin, Zh. Eksp. Teor. Fiz. 58, 1140 (1970) [Sov. Phys. JETP 31, 936 (1970)]. 
[3] A.P. Volodin, M.S. Khaikin, and V.S. Edelman, Pis'ma Zh. Eksp. Teor. Fiz. 26, 707 

(1977) [JETP Lett. 26, 543 (1977)]; U. Albrecht and P. Leiderer, Europhys. Lett. 3, 

705 (1987). 

[4] V.B. Shikin, JETP Lett. 27, 39 (1978); M.M. Salomaa and G.A. Williams, Phys. Rev. 

Lett. 47, 1730 (1981). 
[5] M.M. Salomaa and G.A. Williams, Phys. Rev. Lett. 47, 1730 (1981). 
[6] V.V. Rotkin, and R.A. Suris, Sov.- Solid State Physics 36, 1899 (1994). 
[7] R.D. Averitt, D. Sarkar, and N.J. Halas, Phys. Rev. Lett. 78, 4217 (1997); E. Prodan 

and P. Nordlander, Chem. Phys. Lett. 352, 140 (2002). 
[8] J. Tempere, l.F. Silvera, and J.T. Devreese, Phys. Rev. Lett. 87, 275301 (2001). 
[9] L.P. Gor'kov and D.M. Chernikova, Pis'ma Zh. Eksp. Teor. Fiz. 18, 119 (1973) [JETP 

Lett. 18, 68 (1973)]. 

[10] J. Tempere, LF. Silvera, and J.T. Devreese, Phys. Rev. B 67, 035402 (2003). 

[11] V.B. Shikin and Yu.P. Monarkha, Zh. Eksp. Teor. Fiz. 65, 751 (1973) [Sov. Phys. JETP 

38, 373 (1973)]. 
[12] H. Frohhch, Advances in Physics 3, 325 (1954). 

[13] S.A. Jackson and P.M. Platzman, Phys Rev. B 24, 499 (1981); ibid. 25, 4886 (1982). 

[14] S.A. Jackson and F.M. Peeters, Phys. Rev. B 30, 4196 (1984). 

[15] S. Fratini and P. Quemerais, Eur. Phys. J. B 14, 99 (2000); ibid. 29, 41 (2002). 

[16] F. Lindemann, Z. Phys. 11, 609 (1910); CM. Care and N.H. March, Adv. Phys. 24, 

101 (1975). 

[17] S.N. Klimin, V.M. Fomin, J. Tempere, LF. Silvera, and J. T. Devreese, accepted for 

publication in Solid State Comm. 
[18] H. Ikezi, R.W. Gianetta and P.M. Platzman, Phys. Rev. B 25, 4488 (1982). 
[19] G.E. Marques and N. Studart, Phys. Rev. B 39, 4133 (1989). 

[20] R.P. Feynman, Phys. Rev. 97, 660 (1955); R.P. Feynman and A.R. Hibbs, Quantum 

Mechanics and Path Integrals (Mc Graw-Hill, New York, 1965). 
[21] J.T. Devreese, S.N. Khmin, V.M. Fomin, F. Brosens, Solid State Commun. 114, 305 

(2000). 

[22] C.C. Grimes and G. Adams, Phys. Rev. Lett. 42, 795 (1979). 

[23] V.M. Bedanov and F.M. Peeters, Phys. Rev. B 49, 2667 (1994). 

[24] V.M. Bedanov, G.V. Gadiyak, and Y.E. Lozovik, Phys. Lett. A 109, 289 (1985). 

[25] A. Perez-Garrido, M.J.W. Dodgson, and M.A. Moore, Phys. Rev. B 56, 3640 (1997). 

[26] J. J. Thomson, Philos. Mag. 7, 237 (1904). 

[27] P. Lenz and D.R. Nelson, Phys. Rev. Lett. 87, 125703 (2001). 

[28] D.R. Nelson and B.l. Halperin, Phys. Rev. B 19, 2457 (1979). 

[29] l.F. Silvera, BuU. Am. Phys. Soc. 46, 1016 (2001). 

[30] D.S. Fisher, B.L Halperin, and P.M. Platzman, Phys. Rev. Lett. 42, 798 (1979). 
[31] G. Deville et al, Phys. Rev. Lett. 53, 588 (1984). 



17 



TABLES 



TABLE I. For typical multielectron bubbles, the values of the several physical quantities are 

given at zero applied external pressure. The number of electrons in the bubble (N), the bub- 
ble radius (Rb), the average interelectron distance (d), the surface density of electrons (n^), the 
pressing field generated by the electrons |E| (see Eq. 9), the characteristic energy scale of the 
electron-ripplon interaction [(e|E|)^/(T, cf. formula (16)], and the characteristic frequency of the 
lattice potental {toiat) are given. Compare these quantities to, for example, the maximum density 
(Ri 2 X 10^ cm~^) and the maximum pressing field (?« 3 kV/cm) achievable on a flat helium surface 
over bulk. 



N 




c?(nm) 


ns(cm ^) 


10^ 


0.228 


25.57 


1.529x10^1 


10^ 


4.937 


55.34 


3.265x10^° 


10^ 


106.4 


119.3 


7.025x10" 


N 


|E|(kV/cm) 


e2|E|Va (meV) 




10^ 


138.3 


85.16 


3.891 


10^ 


63.80 


3.884 


1.222 


10^ 


6.350 


0.180 


0.386 
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FIGURES 



FIG. 1. The variational parameter a describing the width of the electron wave function in the 
strong-coupling approach as a function of number of electrons and pressure in the multielectron 
bubble, d is the interelectron separation. 

FIG. 2. For a bubble with N = 10^ electrons, at different pressures the shape of the dimpled 
surface is shown. Electrons are present on the surface, separated from eachother by the lattice pa- 
rameter d. Underneath each electron there is an individual dimple, induced by the electron-ripplon 
interaction. As the pressure is increased, the bubble radius decreases, and the electron-ripplon 
interaction becomes stronger, resulting in a stronger dimpling effect. 

FIG. 3. The phase diagram for the spherical 2D layer of electrons in the MEB. Above a critical 
pressure, a ripplopolaron solid (a Wigner lattice of electrons with dimples in the helium surface 
underneath them) is formed. Below the critical pressure, the ripplopolaron solid melts into an 
electron liquid through dissociation of ripplopolarons. 

FIG. 4. The phase diagram shown in Fig. 3 is extended to reveal the relation of the ripplopo- 
laron Wigner lattice to the Wigner lattice of electrons. These arc distinct, not only in melting 
properties (the ripplopolaron Wigner lattice melts through dissociation of ripplopolarons), but 
also in their location on the phase diagram. The region for the Wigner lattice of electrons without 
dimples -in agreement with the observation of Grimes and Adams [22]- starts at large N and is 
quantum molten by pressurizing the bubble. 
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